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考虑 量化 不 匹配 问题 的 信息 物理 系统 的 安全 控制 ， 
贾 欣 婷 *， 郑 相 超 **， 裴 斌 ? 
(南京 信息 工程 大 学 a. 自动 化 学 院 ; b. 江苏 省 大 气 环境 与 装备 技术 协同 创新 中 心 ， 南京 210044) 
摘 要 :针对 一 类 包含 执行 器 攻击 、 外 部 干扰 和 编码 器 /解码 器 不 匹配 的 信息 物理 系统 (cyber-physical systems, CPS)， 
设计 一 种 新 型 便 棒 控制 器 确保 其 安全 稳定 运行 。 控 制 器 线性 部 分 的 增益 由 优化 H2 性 能 的 线性 给 阵 不 等 式 给 出 ; 控 


制 器 的 非 线性 结构 包含 两 部 分 ,一 部 分 用 于 消除 
估计 来 实现 对 攻击 的 补偿 。 最 后 ,仿真 算 例 乡 
有 效 性 。 

关键 词 : 信息 物理 系统 ; 执行 器 攻击 ; 编码 器 /解码 器 失 配 ; 一 臻 最终 
中 图 分 类 号 : TP273.3 doi: 10.19734/j.issn.1001-3695.2018.11.0890 


余 外 部 和 干扰、 量化 误差 
吉 果 表明 控制 


， 男 一 部 分 利用 对 未 知 执行 器 攻击 参数 的 自 适 应 
器 能 保证 闭 环 系统 的 一 致 最 终 有 界 性 ， 说 明了 所 提 方 法 的 
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Security control of cyber-physical systems considering quantization mismatch 


Jia Xinting’, Zheng Bochaoa pb, Pei Bins 
(a. School of Automation, b. Jiangsu Collaborative Innovation Center on Atmospheric Environment & Equipment 
Technology, Nanjing University of Information Science & Technology, Nanjing 210044, China) 


Abstract: This paper proposed a novel robust controller to ensure the security and stability for a class of cyber-physical 
systems subject to the effects including actuator attacks, external disturbances, and encoder/decoder mismatches. An 
optimized H2 performance in the form of linear matrix inequality solves the linear gain of the constructed controller. The 
nonlinear structure of the control includes two parts: one aims to eliminate the influence of external disturbance and 
quantization error, and the other part compensates the actuator attack by means of adaptive estimations of unknown actuator 
attack parameters. Finally, the simulation results show that the controller guarantees uniform ultimate boundedness of the 
closed-loop dynamical system, which shows the effectiveness of the proposed method. 
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i 可 以 看 做 是 编码 和 解码 的 过 程 ， 编 码 器 和 解码 器 的 量化 参数 
必须 是 相等 ，Yun 等 人 (2 研究 了 具有 输入 量化 和 外 部 扰动 的 

随 着 信息 技术 的 飞速 发 展 及 其 数据 处 理 能 力 的 不 断 提升 ， 不 确定 线性 系统 的 厂 状态 反馈 控制 器 的 设计 。 由 于 硬件 执行 

使 通信 和 控制 系统 越 来 越 紧密 地 联合 在 一 起 ， 信 息 物理 系统 ”的 不 理想 ， 在 实践 中 可 能 会 导致 量化 过 程 中 编码 和 解码 的 量 

GG。 (cyberphysical systems,CPS) 作 为 一 种 新 型 智能 系统 应 运 而 ”化 参数 不 一 致 的 问题 。 针 对 此 缺陷 ，Zheng 等 人 63 研究 了 具 

爷 。。 生 om。CpPS 集成 了 物理 过 程 ， 计 算 资 源 和 通信 功能 ， 本 质 上 有 量化 参数 不 匹配 的 线性 系统 的 也 控制 器 的 设计 。 

三 。 是 具有 控制 属性 的 网 络 , 将 物理 设备 连 上 互联 网 ,但 CPS 更 然而 ,在 CPS 的 分 析 与 控制 的 研究 中 ， 极 少 有 考虑 含 信 
强调 循环 反馈 。 尽 管 CPS 具有 良好 的 广泛 的 应 用 前 景 思 ， 但 “号 量化 ， 特 别 是 量化 编 解码 不 一 致 的 成 果 发 表 。 为 此 ， 本 文 
仍然 有 许多 不 足 之 处 需要 解决 ， 尤 其 是 安全 方面 的 问题 -9。 研究 执行 器 攻击 模式 下 的 含 模型 不 确定 1 的 CPS， 并 同时 考 

于 CPS 使 用 开放 的 计算 和 通信 平台 架构 , 很 容易 受到 虑 量化 不 匹配 问题 ， 通 过 设计 控制 器 以 确保 系统 的 稳定 性 。 
政 对 攻击 ， 攻 击 者 可 以 访问 传 感 和 驱动 计算 平台 ， 操 纵 系统 。 “控制 器 的 线性 部 分 使 模型 不 确定 和 量化 不 匹配 的 系统 实现 鲁 
测量 数据 和 控制 输入 命令 ， 严 重 影响 系统 的 性 能 和 完整 性 ， 棒 丰 性能， 非 线性 部 分 用 于 消除 外 部 干扰 、 量化 误差 以 及 抑 
因此 CPS 中 的 安全 性 比 一 般 计 算 系统 更 重要 中 ,对 于 CPS 的 “， 制 或 抵消 状态 相关 的 执行 器 攻击 对 系统 的 影响 ， 从 而 保证 闭 
可 靠 控 制 问题 ， 针 对 考虑 传感器 攻击 的 CPS，Yucelen 等 人 四 环 系 统 的 一 臻 最终 有 界 性 。 
提出 了 自 适应 控制 器 以 确保 系统 的 稳定 性 。 针 对 考虑 执行 器 i 
攻击 和 噪声 的 _CPS，Xie 等 人 四 提出 了 新 颖 的 切换 观测 器 以 1 预备 知识 及 问题 描述 
确保 系统 的 稳定 性 。 针 对 同时 考虑 传感器 和 执行 器 攻击 的 。 1.1 预备 知识 
CPS，Jin 等 人 (9 提出 了 自 适应 控制 器 以 确保 系统 的 稳定 性 首先 介绍 本 文 会 用 到 的 投影 Paj:RxR > 惨 。 给 出 一 个 
于 数字 通信 信道 在 现代 工业 控制 系统 中 的 广泛 应 用 ， 人 (s+1)070— 0 二 
而 通信 网 络 又 往往 受到 通信 带宽 、 有 限 数据 率 等 困扰 ， 所 以 ”连续 可 微 的 让 函数 为 OO- ai 本 ”， 其 中 0m eR 是 
通常 需要 对 所 测 得 的 信息 进行 量化 处 理 LD。 信 号 量化 的 过 程 。 9eg' 的 投影 范 数 界 ， 忆 >0 是 投影 容 差 界限 ， 对 于 yeR*"， 向 
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量 的 投影 定义 为 码 器 一 侧 获得 量化 测量 信息 qt)/4.(D)) ， 在 传输 编码 过 的 控 
Proj(0,y) = 制 信号 到 解码 器 处 的 网 络 通信 链 路 中 ， 受 到 网 络 攻击 者 恶意 
y if $(0) <0 注入 的 执行 器 攻击 5(1,x()) ， 在 传输 到 执行 器 之 前 的 控制 信 
y if (0)>0,0"(0)y<0 (1) ”号 需要 先进 行 解码 ， 经 解码 器 获得 遭受 执行 器 攻击 的 量化 输 

六 Sm MO) if $0)>0,8"(0)y>0 入 8(aD)， 量 化 反馈 的 信息 物理 系统 如 图 1 所 示 。 


在 向 量 的 投影 基础 上 上， 定义 矩阵 的 投影 


Proj, : Rom x Room 一 > RY™” 


Proj,(9, 了 ) = 
(Proj(col (9),col (了 ))……,Proj(col, (©), col,, (Y))) O) 
中 : Be Rom ， Ye Rw ，eol() 表示 矩阵 4 的 第 ; 列 。 
文中 表示 实数 集合 ，R" 表示 7 维 实数 列 向 量 ， Re”" 表 
示 nxm 维 的 实 和 矩阵 ,47 表示 甜 阵 4 的 转 置 ，4- 表示 矩阵 4 的 


逆 矩 阵 ，X>0CX>0) 表示 X 是 正定 〈 半 正定 ) 和 矩阵， 刷 , 表示 


向 量 的 ? - 范 数 , |, = the tt，p>1, 当 P=% 


时 ， Il =maxeelzs| ， 特 别 的 ， 川 表示 2- 范 数 ， 必 表示 
Frobenius 矩阵 范 数 ，F(z) 表 示 关 的 上 界 〈 下 界 )， 对 称 和 矩阵 
的 对 称 位 置 用 “*” 表 示 。 
1.2 问题 描述 

针对 遭受 执行 器 攻击 的 不 确定 信息 物理 系统 ， 建 立 如 式 
(3) 所 示 的 动态 方程 模型 。 
X(D) =(A+ AAC) xX) + Biilt) + Bo(t) (3) 
其 中 : x(D eR" 表示 系统 的 状态 ，o(D) eR" 是 外 部 干扰 ，A 和 B 
为 已 知 的 具有 适当 维 数 的 常 值 矩 阵 ，A4(7) 是 适当 维 数 的 不 确 
定 参数 矩阵 ，xDsR" 为 遭受 执行 器 攻击 的 控制 输入 ,其 数学 


A(D) =u(D) +6, (t,x() (4) 
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其 中 : .xD):RxR 一 R" 为 攻击 模式 , 是 发 生 在 控制 器 与 执行 
器 之 间 通 信和 链 路 上 的 攻击 。 同 文献 [10] 一 样 ， 本 文 假设 执行 
器 攻击 参数 化 为 6G,x00)=W7TQ9(x(D)) ， 其 中 ，W(n)e Rr 是 未 


知 的 时 变 加 权 和 矩阵 


构 已 知 的 关于 xD) 的 非 线 性 函数 。 从 表达 式 (2) 可 见 ， 如 果 
6.(,") 非 零 ， 则 控制 输入 信号 xD 会 被 错误 〈 或 恶意 ) 的 信号 
2 取代 ， 即 ， 系 统 遭 受 执行 器 攻击 。 
另 一 方面 ， 控 制 系统 的 执行 过 程 越 来 越 多 的 采用 远程 通 
信 或 数字 方式 来 实现 ， 而 数字 通信 网 络 往往 受到 通信 带宽 、 
有 限 数 据 率 等 困扰 ， 所 以 通常 需要 对 所 测 得 的 信息 进行 量化 
处 理 。 量化 器 20O 由 向 最 接近 的 整数 合 入 函数 40 定义 ， 本文 
采用 如 下 的 量化 器 形式 : 
ult 
ou -wosl 2 ] (5) 
其 中 ，j(D 和 人 GD) 分 别 为 编码 侧 和 解码 侧 的 量化 灵 
令 r0D=OO/ 人 GD ,理想 情况 下 ww() 和 从 GD 是 相等 的 , 即 ~=1 。 
然而 在 实际 控制 工程 中 ， 由 于 硬件 执行 不 理想 ， 该 


mwl < 丈 ,| 了 Oo WW ，p(x(D)eR? 是 结 


的 情况 ， r(D) el(r,7), 其 中 ， r+ 和 F 都 是 正 参 数 并 且 7> 
所 以 r=1 也 是 该 情况 下 的 一 种 特殊 情况 。 
综合 上 述 式 (3) ~ (5)， 系 统 (3) 变 为 式 (6) 的 形式 


X(t) =(A+AAG) xX) + BQ)) + Bot) (6) 
其 中 CGO)) 是 受到 执行 器 攻击 的 量化 输入 
= CD 
QD)) = ol a2 2D) 6, (x0)| (7) 


在 信息 物理 系统 (6) 的 运行 过 程 中 ， 控 制 输入 wo) 经 编 


被 控 对 象 


CD) 


解码 器 


图 1 信息 物理 系统 结构 图 
Fig.1 Structure diagram of cyber-physical system 
首先 给 出 以 下 3 个 假设 : 
假设 1 系统 (4,5 ) 是 可 控 的 。 
假设 203 时 变 的 不 确定 矩阵 A4Q) 满足 
AA(D) = FADE ,ADA <7 (8) 
其 中 : 五 和 互 是 已 知 的 具有 适当 维 数 的 常 值 矩阵 ， A() 是 未 
知 的 时 变 矩 阵 。 
假设 303] 存在 已 知 的 参数 s 使 得 外 部 干扰 o(0) 满足 
lo), <e, (9) 
在 给 出 本 文 的 控制 器 设计 之 前 , 首先 引入 如 下 的 五 性 能 
指标 的 概念 [9: 
J (1) = [ {XT (s)Qx(s) +xT (s)KTRKx(s) ds (10) 
此 中 ， 0 和 有 是 给 定 的 对 称 正定 加 权 和 矩阵 ，K 是 控制 器 中 给 
出 的 线性 反馈 增益 矩阵 。 
本 文 用 到 如 下 三 个 引 理 。 
引 理 102 对 于 任意 两 个 4 维 实 向 量 4 与 8， 当 p>1， 
49>1 且 p"'+q7=1 时 ， 下 列 不 等 式 成 立 


le’Blslal, ls 


TS 


qa 


引 理 203 对 于 任意 给 定 适 当 维 数 的 实 和 矩阵 IE 、X 和 7， 
其 中 二 是 对 称 的 ， 当 和 矩阵 FQ 满足 严 (OEO<s7 时 ， 则 
II+XF(DY+YTF() XT <0 
且 仅 当 存 在 一 个 常数 ae>0 使 得 下 列 不 等 式 成 立 
II+eXX +e "YY <0 
引 理 309 根据 矩阵 投影 Proj 的 定义 ， 对 于 给 定 的 
GO' eR"" ， 下 列 不 等 式 成 立 : 

tr(Q—O*) (Proj,,(O, N)—N))<0 


2 ”量化 不 匹配 CPS 的 安全 控制 


针对 上 述 信 息 物理 系统 式 (6) 设计 一 种 新 型 鲁 棒 控制 器 
保证 闭环 系统 的 一 致 最 终 有 界 性 。 其 线性 部 分 的 增益 由 优化 
五, 性 能 的 线性 矩阵 不 等 式 给 出 ; 其 非 线 性 部 分 用 于 消除 外 部 
干扰 、 量 化 误差 以 及 抑制 或 抵消 状态 相关 的 执行 器 攻击 对 系 
统 的 影响 。 构 造 如 下 结构 的 控制 器 

u(D) = Kx(D) +u, (D+v(D) (11) 


Ik 
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其 中 : (A+FRAIB)P+P(A+RHAE)+rKTBP+rPBK < 
a ATP+PA+MIE!E +mPRF P+rPBK +rK'B'P 
Ww) = So ign PB) (12) 
所 以 
A (13) Vx (ATP+PA+IE!E +mPR FP+rPBK +rK’BTP)x 
v(D) = -HW p(x(t)) 

如 果 不 等 式 


ATP+ PA+N'ETE +mPF'FP 


Y GO 满足 自 适 应 律 +1K7BTIP+FPBKE+TO+KTIRK<0 (18) 
WO) = .Proj, (WD), p(x(1)x" PB) (14) 成 并 ,可 以 得 到 , 
VD) <—x’ (PO+ KTRK)x() < fj) 
参数 7 为 正 的 调节 增益 ， Proj, 为 矩阵 的 投影 。 由 于 VC9=V(9=0 以 及 VDsj() ， 可 以 得 到 
Kx(?) 为 线性 部 分 ,用 于 使 不 确定 性 系统 实现 如 式 (10) 的 J(D<V(x(D)<V(x(0)) ， 从 上 述 结论 可 以 看 出 性 能 的 上 界 依 
HH 性 能 。 赖 于 系统 的 初始 状态 x(0) ， 而 在 实际 应 用 中 ， 很 难 精确 确定 
wlD) 与 Y0D) 为 非 线 性 部 分 ,分 别 用 来 处 理 量化 误差 和 外 部 系统 的 初始 状态 , 为 了 克服 这 一 困难 , 可 以 假定 初始 状态 x(0) 
干扰 、 抑 制 或 抵消 状态 相关 的 执行 器 攻击 的 影响 。 是 一 个 满足 E{x(0)x7(0)}=7 的 零 均值 随机 变量 。 通 过 考虑 性 能 
Rs a u i 指标 的 期 望 值 ， 得 到 
为 了 证 明 方 便 ， 定 义 太 = 多 -从 ， “=p )-u, 容易 得 Ee 
到 : 根据 Schur 补 引 理 023， 定 理 1 中 的 约束 条 件 式 (16) 等 
, ), u 价 于 M>P>0, 因此 ，M 的 迹 最 小 化 将 保证 P 的 迹 的 最 小 化 ， 
0 即 系统 素性 能 上 界 的 最 小 化 。 在 式 (18) 两 边 同 乘 以 P-， 
=(A+ RAP)AD + rBu te + HW (ox D+ BoW) 令 F= 户 以 及 K=RKP1 ,根据 Schur 补 引 理 , 可 以 得 到 如 式 (15) 
=(A+ FAE)x(D) +rBLKx() +u, + es + WT (tp(x())] 的 不 等 式 。 因 此 , 通过 最 小 化 正定 和 矩阵 5、M 、 常 规矩 阵 K 、 
+ BO(t) 正 标量 ,可 以 实现 最 小 名 性 能 上 界 , 由 此 定理 1 证 明 完 毕 。 
定理 2 对 于 满足 假设 1~3 的 信息 物理 系统 式 (6), 对 于 


和 WO) =WOD nr Proj, (WO), p(x(t)) x PB i 
WA Mi 所 有 (G0),W(0)e RrxRem ， 由 式 (11) 给 出 的 满足 式 (12) ~ 


定理 1 对 于 满足 假设 1~3 的 信息 物理 系统 式 (6) 并 且 。”“〈14) 的 控制 器 使 闭环 系统 是 一 致 有 界 的 ， 且 xD 入 (最 
5,(4x()=0 ， 当 re(x,7) 时， 存在 对 称 正定 矩阵 5 和 MM， 常规 终 有 界 。 
和 矩阵， 正 标量 mn， 使 以 下 的 优化 问题 有 解 。 


1 Padr'ia, WW+W ) 
min Trace(M) Ix I 元 Ey Cre ) 十 ?7 (多 十 站 (19) 
和 [Wo < tas PYdr'ds + (W + We 2 (20) 
= 
* OO- 0 <0 (15) 2 
上 4。 KR 中 dA (QtKIRK) ,d= Wt 
M I so 和 a 
1 5|>0 (16) 2771(W+ 记 ,DW ， 六 ,eR ，W eR 是 投影 的 范 数 界 。 
其 中 : T=AP+PA7+mFF7+rBK+rK7B7 。 当 k=RP21 时 ， 则 由 式 证 明 考虑 李 雅 普 诺 夫 函数 V=x(0)7Px(D)+1tr(W7(DW0)) ， 


可 


(11) 给 出 的 满足 式 (12) 的 控制 器 使 闭环 系统 x0 轨迹 渐 前面 的 证 明 类 似 定理 1，Y) 沿 系统 轨迹 的 时 间 导 
近 收 敛 于 原点 并 满足 了 性 能 ， 可 以 通过 最 小 化 矩阵 W 的 迹 
来 达到 最 小 的 丈 性 能 上 界 。 

证 明 ”考虑 李 雅 普 诺 夫 函数 V=xOD7Pxt) ，VYG) 沿 系统 轨 
迹 的 时 间 导 数 为 


数 是 


V=iPrtx Prt 2mir(WrW) < 
—x' (QO+KTRK)x+2rxrPB(u WT o(x(t))) 
—27 01r[WTProj, (W, p(x(D)) x PB)] + 277 1tr (WTW) 


V=X Px+x’ Pi 根据 引 理 3 可 得 
=xi'[(A+ FAIE, )P+P(A+ 五 人 已 )+FK7TBTIP+FPBK]x 2rxT PB(1 WT 9(x())) 一 27 tr[WT Pr oj ( 帮 ,o(x(D))xrPBD)] Es 
+2rx’ PB(u, +e, )+2x7PBmw(D) 2rytrWro(x(O) x PB] 
根据 引 理 1， 假设 2、3， 因 为 不 等 式 =2(r—r +7) AW Proj, (W, p(x(D)x’ PB)] = 
u us uu 2rpetr[(W —WY (Proj, W, p(x(D)x’ PB) — p(x()xTPB)] 
le =aq — |-ul = (4(—)-——) sn. 12 本 本 和 
He 四 He KH ls +2(r —7T)Wtr[WT Pro (W, p(x(1))xT PB)] < 
可 得 2(r —7)utr WT Pr oj, (W, p(x(t) x PB)] 
2rx7 PB(u, +e,)+2x’PBo < 所 以 
2rxT PBu, + 2r|x’PB +2|x’PB| lol, < 
1x u rlx | lex | llx [i lowl (17) Vz 之 = 这 (O+ KTRK)x+ 2(r— 2 WiW) 十 277- py WW) Ee 
2rx’ PBu, + 2rl|x’PB| 全 +2|x7PBl sw -a x 半 苞 
将 式 〈12) 代入 式 (17) 可 得 寻 此 ， 在 Q 集合 之 外 ，V<0。 
2rx7 PBu, +2r|x" PB), 入 +2|xzrPBl 2, <0 Q={(xD, WD) eR xRe" :OFS 3, WO < 2} 


根据 引 理 2 可 得 


其 中 9 = Ya/dt ， 品 = 态 + 刻 。， 这 就 证 明了 定理 2 闭环 系统 的 
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致 有 界 性 现 令 人 满意 的 系统 性 能 。 
因为 30 
| 
hn PDF +17 Wl < oo (21) 20h 和 
| | 3 
其中 :wm = (P+ ， 式 (21) 可 得 lz < Se i _ jo 
WW <woww ， 所 以 x() 和 及 CD 是 最 终 有 界 的 ， 由 此 定理 2 证 明 oF 
完毕 。 ol 
3 。 实验 结果 与 分 析 JR 
ee ee A “0 0.2 0.4 0.6 0.8 1 
下 面 给 出 一 个 仿真 实例 来 验证 本 文 所 提 方 法 的 有 效 性 。 dings 
考虑 如 式 (6) 形 式 的 飞行 器 动态 系统 的 线性 化 模型 0， 其 中 ， 图 2 系统 状态 响应 曲线 图 
系统 状态 xD =[B(,p(D,YODT 包含 侧 滑 角 ( deg )、 滚 转 率 ( deg/s ) Fig. 2 The response curves of system states 


以 及 偏 航 率 ( deg/s )， 控 制 输入 wD =[Guz(D,5u 包含 
遍 转 ( deg )， 系 统 的 矩阵 参数 如 下 : 


-0.025 0.104 -0.994 0.122 -0.276 
4=| 574. ,B=| -53.61 33.25 
195.5 


转 ( deg ) 以 及 方向 舵 1 


574.7 0 
16.20 0 


0 
0 


n= 0 可， 


十 


仿真 中 考虑 三 种 情 
a) 为 了 说 明 量化 不 


机 a 


匹配 问 


副 翼 侦 


一 529.4 


五 =0 1 1 


题 对 系统 稳定 性 


行 和 
和 矩阵 
制 器 


对 


HH 


制 器 


QO=0.01L, ， 


P= 


Simulink 对 


和 


图 


[12] 的 控制 


系统 


参数 


偿 器 wD) 的 
言 息 物 理 系统 进行 仿真 。 
c) 针 对 上 述 修 改 为 不 匹配 的 量化 灵敏 度 参数 ， 
里 1 中 的 Zr 求解 给 出 对 称 正定 矩阵 P 和 控制 器 增益 K， 控 
类 器 wl?) 如 (12) 的 形式 进行 仿真 。 

根据 文献 [12] 中 的 定型 
R=0.11;， 求 得 
2.8242 0.0971 ey 


化 灵敏 度 参数 匹配 时 的 信 
文献 [ 
中 的 非 线性 补偿 器 w() 采 | 
b) 同 上 述 的 对 称 正定 矩阵 了 、 控 制 器 增益 K 选 
将 量化 灵敏 度 参数 修改 为 不 


P 和 控制 器 增益 


EB 式 相 同 ， 


息 物 理 系 统 的 仿 


E 的 影响 ， 先 进 


， 对 称 正 定 


中 的 补 
情况 1 


=| 0.0971 
-0.0102 0.0001 


0.0043 ”0.0001 | ， 
0.0006 


12] 定 理 


选取 量化 灵敏 度 参 


6, (t,x(t)) =[0.5sint cos(Q2DT 0.2sin(B(1)) cos(p(t)) 。 


数 Ha 


= 从 =0.00 ， 


EE 3， 选 取 参 数 矩 


天 二 35.4708 1.4384 -0.0406 | 
| 6.4268 0.4774 0.2267 | 


文献 [12] 的 形式 。 


取 以 及 补 
匹配 时 ， 


由 本 文 定 


阵 


假设 执行 器 攻击 为 
运用 Matlab / 


系统 (6) 进行 仿 


， 得 到 系统 的 响 


应 


线 如 图 2 


3 所 示 。 


趋 近 于 稳定 。 
情况 2 


发 生 不 


匹配 现象 ， 


可 以 看 出 ， 
器 使 信息 物理 


当量 
系统 


选取 量化 灵敏 度 参 
运用 Matlab/Simulink 进行 仿真 ， 得 到 系统 的 响应 
和 5 所 示 ， 念 真 结果 看 


闭 i 


化 器 灵敏 度 参 数 


的 状态 和 控制 输入 收敛 


> 数 分 别 为 Ha= 


匹配 时 ， 文 献 
于 原点 ， 


=0.06 和 =0.1， 


线 如 图 4 


上 ， 由 于 编码 器 和 解码 器 量化 灵敏 度 


稳定 现象 ， 文 献 [12] 的 方法 失效 。 


ro A 
所 示 ， 仿 真 结 果 表 明 设 计 的 控 


情况 3 根据 本 文 定型 
0.9312 0.0369 0.0030 
oe 0.0024 0.0008 | ， 
0.0030 0.0008 0.0007 


pie 


运 


1, 


本 文 定理 2 中 的 控制 


可 求 得 


he -0.3813 
天 Es 


| 9.5125 5.1038 


的 响应 


[于 
x 


环 系统 状态 是 发 散 的 ， 系 统 呈 现 不 


一 1.3382 
4.9255 


器 设计 方法 ， 并 经 MATLAB 
线 如 


多 


6 和 7 


制 器 在 信息 物理 


系统 发 生 量化 


灵敏 


度 参 数 不 


匹配 、 尹 


部 干扰 


以 及 模型 不 确定 性 


时 ， 仍 能 实 


u(t) 


x(t) 


u(t) 


40 


305 


20. 


区 | 


0.2 


3 


04 0. 6 
time/s 


系统 控制 输入 响应 


1 线 


图 


Fig.3 The response curves of control inputs 
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图 4 系统 状态 响应 曲线 图 


Fig.4 The response curves of system states 
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x 10° 


区 | 
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time/s 


系统 控制 输入 响应 


线 


图 


Fig.5 The response curves of control inputs 
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图 6 系统 状态 响应 曲线 图 


Fig.6 The response curves of system states 


30F 


25 
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time/s 


7 系统 控制 输入 响应 曲线 图 

Fig.7 The response curves of control inputs 
从 上 面 三 种 情况 的 仿真 效果 对 比 表 明 ， 文 献 [12] 设 计 的 
1 器 能 确保 在 发 生 情况 1 的 信息 物理 系统 ， 即 量化 匹配 的 
息 物 理 系 统 ， 实 现 稳 定 ， 但 当 发 生 量化 不 匹配 时 ， 由 仿真 
4 和 5 可 以 看 出 文献 [12] 的 控制 器 已 经 无 法 确保 信息 物理 
系统 的 稳定 ， 严 重 时 甚至 导致 系统 失 稳 。 由 情况 @ 的 仿真 可 
以 看 出 , 本文 设计 的 控制 器 能 有 效 克 服 量 化 不 匹配 等 的 影响 ， 
确保 信息 物理 系统 的 稳定 ， 充 分 验证 了 本 文 控制 算法 的 有 效 
生 和 优越 性 。 


4 ”结束 语 


本 文 研究 了 存在 执行 器 攻击 、 外 部 干扰 和 量化 不 匹配 的 
信息 物理 系统 的 安全 控制 问题 。 设 计 的 控制 器 由 线性 和 非 线 
性 两 部 分 组 成 ， 甚 中， 线性 部 分 用 于 使 模型 不 确定 性 系统 实 
现 五 性 能 ; 其 非 线 性 部 分 用 于 消除 外 部 干扰 、 量 化 误差 以 及 
抑制 或 抵消 状态 相关 的 执行 器 攻击 对 系统 的 影响 ， 以 确保 闭 
环 系统 的 一 致 最终 有 界 性 。 仿 真 结果 体现 了 文中 算法 的 有 效 
性 与 优越 性 。 未 来 的 研究 工作 中 ， 将 重点 解决 发 生 通信 中 断 
和 时 延 等 情形 下 的 信息 物理 系统 的 安全 可 靠 控 制 问题 。 
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